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NOTATION 


b Angular orientation at time t, not necessarily along the most direct rota- 

tional path; also used to indicate the body system 

bo Initial angular orientation 

bf Final angular orientation 

b s Angular orientation, for following the most direct rotational path 

D u _> v 3x3 element direction cosine matrix (DCM), for transforming a coordinate 

of the u-system into the v-system. 

e Exponent of disturbance function 

£ - {e x ,e y ,e z } T Unity vector, indicating the direction of the Euler axis for rotation from an 

initial system to a final one, expressed in the initial system 

M u _>v 4x4 element quaternion matrix of the quaternion for rotation from u to v 

* 

M u _^ v Transmutated quaternion matrix 

q u _» v = ( < Il» c l2*Q3’ c l4} T Quaternion for rotation from u to v, uniquely describing the orientation of 

v into u 

T Transpose of a matrix; also used to indicate a column vector 

t Time 


tq Initial time 

tf Final time 

£ b = {x b ,y b ,z b } T Orthogonal right-hand object-body coordinate system; in aeronautical 

applications the x b -axis is pointing forward along the aircraft body main 

axis, the y b -axis is pointing outward through the right wing, and the 
z b -axis is pointing downward 

x* = {x\y\z*} T Orthogonal right-hand Earth-based inertial reference system, with the 

x*-axis pointing to the north, the y'-axis to the east, and the z'-axis down- 
ward, toward the center of Earth; in the experiments the x'-axis coincides 
with the direction of the viewing axis 


m 
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*qO Like the initial system bo, fixed with respect to the inertial system, with 

the xq°-axis aligned with the Euler axis and the yq°-axis in the x b0 -o-y b0 
plane 

xql ;x^ 2 ;xq 3 Intermediate systems to obtain the plane in which to rotate £D b in order to 

be disturbed 

“Angular distance” between orientation b and bf, i.e., the angle of rotation 
for the remaining most direct path rotation from b to bf 

Total angle of rotation for the most direct path rotation from initial orien- 
tation bo to final orientation bf 

“Rotational deviation” from the most direct path; it is the amount of rota- 
tion of the quaternion Aq 

Averaged rotational deviation from the most direct path over the interval 
to < t < tf 

Azimuth angle and elevation angle, respectively, which specify the orien- 
tation of £0o in the bo-system 

Quaternion for the rotational deviation, i.e., for rotation from b s to b 

v Angle of disturbance; it is the angle over which £ 0 b is rotated in order to 

be disturbed 

vo Angle of disturbance at time to 

Azimuth angle and elevation angle, respectively, which specify the 
orientation of £oq° in the qo-system 

\|/,0,(|) Euler angles for yaw, pitch, and roll, respectively, specifying the orienta- 

tion of body with respect to the inertial system 


a 

ao 

P 

P 

Y.5 


^ Angle of plane in which ojb is rotated to be disturbed 

©0 = fi>o Vector of angular rotation, for the most direct rotational path between ori- 

entation bo and bf, expressed in the bo-system 

too - Ifflol Constant rotational speed 

to b Disturbed vector of angular rotation; it is the vector of rotation at which 

the object proceeds rotating at time t, expressed in the b-system 


IV 


Vector of angular rotation, for the most direct rotational path between ori- 
entation bo and bf, rotating at reduced angular rate, expressed in the 
bo- system 

Vector of angular rotation, for the remaining most direct rotational path 
between orientation b and bf, expressed in the b-system 

The vector expressed in the qo-system 

4x4 element matrix, containing the components of used for comput- 
ing the rotational path 

4x4 element matrix, containing the components of (2s, used for 
computing the most direct rotational path 




AN ALGORITHM FOR THE SYSTEMATIC DISTURBANCE 
OF OPTIMAL ROTATIONAL SOLUTIONS 

Arthur J. Grunwald 1 and Mary K. Kaiser 
Ames Research Center 


SUMMARY 


An algorithm for introducing a systematic rotational disturbance into an optimal (i.e., single axis) 
rotational trajectory is described. This disturbance introduces a motion vector orthogonal to the 
quaternion-defined optimal rotation axis. By altering the magnitude of this vector, the degree of non- 
optimality can be controlled. The metric properties of the distortion parameter are described, with 
analogies to two-dimensional translational motion. 

This algorithm has been implemented in a motion-control program on a three-dimensional graphic 
workstation. It supports a series of human performance studies on the detectability of rotational trajec- 
tory optimality by naive observers. 


1. INTRODUCTION 


This paper describes an algorithm for generating kinematically suboptimal (“warped”) rotational 
trajectories. First, the basic idea of creating a suboptimal trajectory is demonstrated with a two- 
dimensional translatory analogy. Next, the mathematical formulation for describing the rotation of an 
object is given for three methods: (1) Euler angles, (2) direction cosine matrix, and (3) quaternions, and 
the relation between these methods is discussed. Then, the rotational equivalent of the straight path for 
translatory motion is defined and ways for computing the rotational path and the deviation from the 
direct path are delineated. Finally, the method used for creating a reproducible “smoothly warped” sub- 
optimal trajectory is outlined. This method allows both the shape of the trajectory and the magnitude of 
the deviation from the optimal path to be defined by two independent parameters. A flow chart is pre- 
sented which summarizes the computations performed to create displays employing this algorithm. 


'Tcchnion, Haifa, Israel. 



2. DEVIATION FROM THE OPTIMAL PATH TRAJECTORY FOR 

TRANSLATORY MOTION 


2.1 Method for Creating a Suboptimal Trajectory 

The basic idea of creating a suboptimal trajectory can best be illustrated with a simple, two- 
dimensional translational analogy. Consider two points on a plane as illustrated in figure 1. An object 
moves from the initial location bo, at time to, to the final location bf, at time tf, at a constant velocity, 
Vo. The kinematically optimal trajectory between bo and bf is a straight line; the traveled distance 
along this line is Ro- At time t, the object is at location b, which is not necessarily along this most 
direct path. However, the remaining most direct path from b to bf is again given by a straight line 
which connects these two locations. If the velocity vector, V, points from b to location bf, the object 
would proceed from t until it reaches tf along this remaining most direct path. However, in order to 
produce a remaining suboptimal trajectory, a disturbance is introduced by rotating V over the angle of 
disturbance v, to obtain the disturbed vector, Vj. The angle v is chosen to be a simple exponential 
function of the range R between b and bf, according to 


tan v(t) = 


R(t) 


i R °, 


tan Vq 


( 1 ) 


where the exponent e determines the characteristics of the disturbance and therefore the trajectory 
shape. The parameter vo, which is the angle of disturbance at to, determines the magnitude of the 
disturbance. 

For e > 1, v(t) will reduce to zero quickly, such that the trajectory will be curved mainly at the 
beginning and straight toward the end, while for 0 < e < 1 the opposite is the case. Typical trajectory 
shapes for these two ranges of e are shown in figure 2. 

2.2 Measures of “Suboptimality” of Motion 

Various measures can be considered for specifying the “degree of suboptimality” of motion. The 
first one is the difference between the traveled distance along the nondirect path, and the shortest dis- 
tance, Ro- For constant velocity Vo, this is equivalent to the difference between the actual travel time 
(tf - to) and the most direct path travel time, Ro/Vo- Since this difference will be small relative to the 
total travel distance (or travel time), this measure will be fairly insensitive for trajectories with small 
deviations. 
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A second candidate for a measure of indirectness is the averaged angle of disturbance, computed by 


v = 


1 

(tf - to) 



v(t)dt 


( 2 ) 


This method has the disadvantage that only the disturbance and not the actual deviation from the straight 
path is considered. Thus, for a straight trajectory section (like the final section of a trajectory with 
e > 1), v is zero and does not add to the measure of indirectness, although the deviation from the most 
direct path definitely exists (fig. 2a). 

A third measure of indirectness is the averaged deviation from the straight path. The way in which 
this score is computed is shown in figure 3. It is assumed that the object travels along the nondirect path 
with a constant velocity Vo. At each instance of time t, the distance d(t) between a point on the tra- 
jectory b and the equivalent position on the optimal path b s is computed; b s is the position on the 
optimal path at which the object would arrive at time t when traveling with a reduced speed of 
V s = Ro/(tf- to) < Vo- The positions bandb s are computed by solving the differential vector equations 

Vt> = (^v^o (3) 


with initial conditions 


Xb( t 0) = 3b s ( t 0)=^b 0 


where 


*b = 


*b c 


and the distance d(t) is 


d(t) = lx b (t)-Xb s ( t )l 


(4) 


The averaged deviation from the optimal path over the interval to ^ t < tf is then computed by 

r*f 


d = 


-j_r 


d(t)dt 


(5) 
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3. DEVIATION FROM THE OPTIMAL TRAJECTORY FOR ROTATIONS 


3.1 Definition of Coordinate Systems and Object Angular Rotation 

The description of the angular rotation of an object in space involves rotational transformations 
between coordinate systems. Two basic systems are defined: an Earth-based inertial reference system, 
i, and an object body coordinate system, b. The Earth-based system, defined according to the aeronauti- 
cal convention, is a right-hand system with its x*-axis toward the north, its y‘-axis to the east, and its 
z'-axis pointing downward toward the center of Earth. In creating displays on the workstation, the 
inertial system is chosen to coincide with the observer’s viewing system, with the x>-axis being the 
viewing direction. The object is described in the body system. For aerospace applications, for example, 
this system is attached to the aircraft body, with the x b -axis pointing forward along the aircraft body 
main axis, the y b -axis pointing outward through the right wing, and the z b -axis pointing downward. 
(NOTE: Coordinate systems other than the North-East-Down convention are frequently used by non- 
aerospace disciplines. However, by means of simple geometrical computations, coordinates expressed 
in one system can be transformed to another; the generality of the algorithm presented in this paper is 
not affected.) 

The angular orientation of the body system with respect to the Earth system can be described in a 
variety of ways. The first one is by a set of three Euler angles, which define three successive rotations. 

In aerospace applications these angles are the yaw angle \y, which sets the heading plane; the pitch 
angle 0, which sets the aircraft elevation angle with respect to the horizontal plane; and the roll angle <)>, 
in this order (fig. 4). A coordinate defined in the inertial system, x* = (xky^z'} 7 , can now be expressed 
in the rotated body system, x b s {x b ,y b ,z b } T , by 

x b = [(l)][0][\|/]x i = D{ V ,0,<j)}x i = D i ^ b x i (6) 


where 

c\|/ s\(/ 0 c0 0 — s0 10 0 

[\|/]= — s\|/ c\|/ 0 ; [0] = 0 1 0 ; [<)>]= 0 c<J> s0 (7) 

0 0 1 s0 0 c0 0 — s<> c<j> 

where cands denote sine and cosine, respectively, and D is the nine-element direction cosine matrix 
(DCM). The superscript T denotes transpose, since a column vector is specified. Likewise, by the 
inverse transformation the coordinates of the b-system can be expressed in the i-system by 

x' = [y] T [Of [<|> 1 T x b = D T {y,0,(j) }x b = D b ^x b (8) 
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where the superscript T denotes the transposed matrix and, since D is an orthonormal matrix, trans- 
pose and inverse are identical. 

Thus, the second way of expressing the angular orientation of the b-system in the i-system is by 
means of the nine-element DCM. A third way is through the use of quaternions. Euler’s theorem states 
that, regardless of the initial orientation (defined by the system u) and the final orientation (defined by 
the system v), it is always possible to find one axis about which the object can be rotated to bring it 
from orientation u to v. The orientation of this axis of rotation, expressed in the initial system u, is 
given by the unity vector £ = {e x ,e y ,e z } T , and the amount of rotation by the angle a. The quaternion 
for rotation from u to v, which uniquely describes the orientation of v with respect to u, is given by 

cos a/2 
e x sin a/2 
e y sin a/2 
e z sin a/2 


Since £ is unity, it follows from equation (9) that 

iq! = jq j + q 2 + qj + <4) =1 0°) 

which is a unique property of the quaternion, and 

lal = 2 cos -1 qi = 2 sin -1 ^ + q^ + q^j ; 0<a<180° (11) 

It also follows from equation (9) that the inverse quaternion, for rotation from v to u (i.e., the quaternion 
which describes the orientation of the u-system in the v-system), is given by 


-1 

9v-»u ~ Su->v 


which simply means that the direction of rotation about the Euler axis is inverted. 
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It is very useful to find the quaternion for a sequence of rotations. Suppose that the angular orien- 
tation of system v with respect to u is given by q tl „> v and the orientation of a third system, w, with 
respect to v by q v _* w . Then, the orientation of w with respect to u is given by 

Su— >w = { My— ^ w)Qu— >v = {^u— >v}*5y— >w (13) 


with 


qi 

-q2 

-q3 

”*l4 


qi 

-q2 

-q3 

“^4 

q 2 

qi 

q4 

-q3 

; M* = 

q2 

' qi 

1 

-q4 

q3 

q3 

1 -q4 

qi 

q2 


q3 

1 q4 

qi 

-q2 

q4 

! ^3 

-q2 

qi 


q4 

! -q3 

q2 

qi 


(14) 


where M is the quaternion matrix, composed of the elements of q, and M* is the transmutated quater- 
nion matrix, obtained by transposing the vector kernel (or minor) of the first element, indicated by the 
dotted partition. 


By using equation (13) in a sequence of rotations, the quaternion q^^ for rotation from inertial 
system i to body system b can be found as a function of the Euler angles. Following equation (9) the 
quaternions for the yaw, pitch, and roll rotations are given by 


£i-»u - 


c\|//2 


_ c0/2 " 


c<J)/2 

0 

* S.U— >v 

0 


s<()/2 

0 


S0/2 

0 

_ S V /2 _ 


0 


0 


(15) 


respectively, where u and v denote the intermediate stages. Then, qj^b is given by 

! 9i— >b = ( ^v— > 5 } { M u __ >v }Qi — mi (16) 

and after evaluating equation (16) with equations (14) and (15) 

c(|)/2 c0/2 c\j//2 + s(|)/2 s0/2 sv|//2 
s<J)/2 c0/2 c\j//2 - c<j)/2 s0/2 si|r/2 

q;_>b = (1 7 ) 

- c<|>/2 s0/2 c\i//2 + s<)>/2 c0/2 sy/2 

— s<j>/2 s0/2 c\j//2 + c(J)/2 c0/2 sy/2 
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The DCM Di_>b can be expressed in terms of as follows: 


l-*i— >b — 


2 . 2 2 2 

q 1 +q 2 -q 3 - £ l4 

^(^2^3 “ qiQ4) 
2(q2Q4 + qiQ3) 


2(q 2 q 3 + qiq4> 

2 2,2 2 

q i - q 2 + q 3 q 4 

2(q4Q3 - qiq2) 


2(q2q4 - qiq3) 
2(q4Q3 + qiq2) 

9 9 9 2 

qi - q 2 - q 3 + q 4 


( 18 ) 


3.2 Rotational Equivalent of a Translatory Straight Path 

It is clear from the definition of the Euler axis that the rotational analogue of a straight path between 
two points in the translatory case is a rotation about the Euler axis. This rotation brings the object from 
an initial orientation to a final one along the most direct path. The total angle of rotation about this axis 
is equivalent to the distance between two points in the translatory case. Equivalently, this means that 
the instantaneous axis of rotation (or vector of angular velocity) is fixed in space. This also means that 
each coordinate of the object will describe a path, which will be the great-arc of a circle, located in a 
plane perpendicular to this fixed axis of rotation. For points of the object located on this axis, the radius 
of the circle will be zero. Mathematically, this means that the coordinates of these points are not 
affected by the rotational transformation. 

In contrast, when the rotation is not proceeding along the optimal trajectory, the instantaneous axis 
of rotation will not be fixed in space and will perform a “wobbling” motion, somewhat like the nutation 
of a gyroscope. 

3.3 Computation of the Rotational Path 

Similar to translatory motion, in which the path is obtained by integrating the instantaneous velocity 
vector, the rotational path is obtained by integrating the components of the vector of angular velocity. 
For the rotational motion of the body system b, this vector is given by 

ffliUbW* {0) x ,(0y,(0 z } T (19) 


where the subscript i-»b indicates the rotation of the b-system with respect to the i-system, and the 
superscript b indicates that this vector is expressed in the rotating b-system. The rotational path is then 
obtained by solving the differential equation 

qi_>b(t) = ( Q(t)}qi->b(0 ; qi->bOo) = q^ (20) 


where bo indicates the orientation of b at time to and 


7 



Q(t) = 1/2 


0 

-co x 

— 0)y 

-to z 

co x 

0 

(0 Z 


0) y 

-®z 

0 

co x 

co z 

CDy 

-co x 

0 


( 21 ) 


3.4 Computation of the “Deviation” from the Optimal Rotational Trajectory 

K k 

For the optimal rotational trajectory, the vector ffi,_^ b (t) = = ffio = constant will be fixed in 

space and will coincide with the Euler axis, coo =l£Qol is the constant rotational speed and ao is the 
total angle of rotation, equivalent to the shortest distance, Ro, in the translatory case. The computation 
of the “deviation” from the optimal trajectory is analogous to the translatory case (fig. 3). At each 
instant of time the quaternion for the orientation b,q i _ >b (t), and the quaternion for the corresponding 
orientation along the optimal trajectory b s ,qi_> b (t), are computed. b s is the orientation on the optimal 
trajectory which would be reached at time t when rotating at a reduced angular rate of co s = ao/(tf- to) 
< coo- The quaternion q[^, bs (t) is obtained by solving the differential equation 

£i-»b s 0) = {^s^£i->b s W » Si — >b s (^0^ = Qi-»b b (22) 


where Q s is fixed and generated according to equation (21) with 

a 0 


co s = 


(tr-t 0 )«o 


®0 


(23) 


The “rotational deviation” from the optimal trajectory (which is equivalent with the deviation d 
from the straight path in the translatory case) is given by the amount of rotation |3 of the quaternion 


-1 


Aq(t) = q bc _>b(0 = {Mi_^ b (t)}q. , (t) 


(24) 


where q^ is computed from Qi^b s by equation (12) and p is computed with equation (1 1). The 
averaged deviation angle over the rotation interval to ^ t < tf is then computed by 


0- 


Of 




(25) 
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3.5 Introduction of the Disturbance 


Consider that the b-system rotates from initial orientation bo at time to to final orientation bf at 
time tf, not necessarily along the optimal trajectory. At time t the orientation of the b-system is b. 
However, analogous to the translational case (in which the direction of the remaining most direct path is 
V}, in the rotational case the remaining most direct path from time t to tf onward is given by the Euler- 
axis rotation from b to bf 


qb->b f (t) = {Mi^ bf }q b ^i(t) = {M^lq^hO) (26) 


where qj_^ b is computed from qj_> b by equation (12), qi_> b is computed by solving the differential 
equation (eq. (20)), and Mj_* bj is given by the final orientation of b in i. For 0 < a < 180°, the orienta- 
tion of the Euler axis for rotation from b to bf (expressed in the b-system) is specified by 


e(t) = 


02/s 

q 3 /s 

q4/s 


s = 


2 2 2 

q 2 + q 3 + q4 


1/2 


= sin(a/2) 


lb-»bf 


(27) 


and the “angular distance” from btobf is a = 2 sin _1 s (eq. (11)). If the instantaneous vector of rota- 
tion ®i^ b (0 = £> b (t) (the subscript i— »b is henceforth omitted for clarity) would be chosen along g(t), 
the rotation would proceed from t onward until it reaches tf along a remaining most direct path. How- 
ever, in order to produce a remaining suboptimal trajectory, a disturbance is introduced to deviate £fi b (t) 
from £(t). As in the translatory case, the deviation of ffl b (t) from £(t) should be a function of the angu- 
lar distance a(t). The vector fO b (t) is deviated from e(t) by rotating it over the angle v(t), which is 
chosen to be a simple exponential function of a(t) according to 


I 


tan v(t) = 


a(t) 


\e 


a 0 


tan Vq 


(28) 


where Vo is the disturbance at time to, and e is an exponent. The parameter vo determines the aver- 
aged deviation from the optimal trajectory and thus the magnitude of the disturbance, and the exponent 
e determines the shape of the disturbance. Since b globally approaches bf, a(t) will decrease mono- 
tonically such that the time-history of v(t) (and thus the shape of the trajectory) will be determined by 
the exponent e. For e > 1, v(t) will reduce to zero quickly such that the trajectory will be “warped” 
mainly at the beginning. For 0 < e < 1 , the greatest warping occurs toward the end of the trajectory. 

It should be noted that, although the angle of disturbance v(t) is defined, the plane in which is 
rotated is still undefined. A sequence of rotational transformations is required to choose a plane which 
ensures a reproducible, smoothly warped rotational path. The method by which this plane is chosen is 
described below. 
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The first rotational transformation that is required to introduce the disturbance involves a trans- 
formation to the qo-system. Like the bo-system, the qo-system is also fixed with respect to the inertial 
system, but it has the x^-axis aligned with the fflo-axis, and the y9 0 -axis in the x^-o-y 150 plane 
(fig. 5). Since the orientation of fflo in the bo-system is given by its azimuth angle y and elevation 
angle 5, the fixed DCM for rotation from bo to qo is obtained by a successive yaw rotation by the 
angle y and a pitch rotation by the angle 5, following equations (6) and (7), according to 



c5 

0 

-s5 

cy 

sy 

0 


c8cy 

c5sy 

-s8 

Db o-*io - 

0 

1 

0 

-sy 

cy 

0 

= 

-sy 

cy 

0 


s5 

0 

c5 

0 

0 

1 


s5cy 

s5sy 

c5 


First, fi2 b (t) = cooe(t) for the remaining most direct path from time t until tf, as computed by equa- 
tion (27), is transformed from the b-system into the qo-system by 


ffl q 0(t) = [D b _ X j 0 lco b (t) 


(30) 


where 




(31) 


The first two matrices on the right-hand side of equation (31) can be precomputed before starting the 
motion-control program, and the third matrix is computed from q b _>j = q^ b by solving equation (20) 
and using equation (18). 


The orientation of £>1° in the qo-system is shown in figure 6. If the rotation were totally along the 
optimal trajectory, 0)9° would coincide with £Oo and thus with the x^-axis. However, its deviation 
from this axis can be expressed in the qo-system by the azimuth angle x and the elevation angle A 
system qi is now defined, with the x 4 ! 1 -axis along &) 4 ! 0 and the y 4 ! 1 -axis in the x^-o-y 4 * 0 horizontal 
plane. Similar to equation (29), the DCM for rotation from qo to qi is given by 


c^c% c^s% -s£, 


D qo-*u “ 


-sx cx 0 
s^cx s^sx c^ 


(32) 


Since 
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m q ° = {co x ,co y ,(o z ) T ; lfij q o| = (o 0 


(33) 


it follows from the geometry of figure 6 that 

cos % = co x /d ; cos ^ = d/co 0 

sin x = 0) y /d ; sin £ = -G) Z /CO0 

where 


C0 z /0)o 

0 

d/to 0 


d = |to^ + o)y 


1/2 


Substituting these expressions in equation (32) yields 


D, 


qo-*h 


co x /co 0 

to y /co 0 

-co y /d 

co x /d 


-CO x G3 z /(dCOo) -G> y co z /(dcoo) 


(34) 


(35) 


(36) 


Equations (35) and (36) show that Dq ) _ Kll can be computed simply from the components of co4 0 in the 
qo-system. 


Next, the q2-system is defined now rotated about the x^-axis by the angle <j) d (fig. 7). This angle 
is randomly chosen from a look-up table, but remains fixed throughout a trajectory calculation. The 
DCM for rotation from qi to q2 is then given by 


D qi-*l2 - 


1 0 0 
0 c()) d s<> d 

0 -s<}> d c<J> d 


(37) 


The plane x<l 2 -o-yq 2 is the required plane in which the vector & is rotated to be disturbed (fig. 7). 
Last, the system q3 is defined, rotated about the z4 2 -axis by the angle of disturbance v(t), and com- 
puted by equation (28) (fig. 7). The DCM for rotation from q2 to q3 is then given by 
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cv sv 0 


Dq 2 -xi3 = - SV CV 0 (38) 

0 0 1 

With the DCMs previously computed with equations (31) and (36)-(38), the DCM for rotation from 
b to q 3 can be computed according to 

The x^-axis is now aligned with the direction of the “disturbed” vector CQjj. This vector, expressed in 
the b-system, is then computed by 

< 40 > 

where 

£2 q 3s[co 0 ,0,0] T (41) 

2.5 Computer Implementation 

The motion-control program involves (1) precomputation of matrices, vectors, and quaternions, 
which remain constant for all created trajectories, and (2) real-time computations for realizing the rota- 
tional path and computation of the deviations from the most direct path. This motion-control program is 
incorporated into an experimental control program for psychophysical studies of human observers sen- 
sitivity for trajectoiy optimality. This experimental control program presents the motion stimuli, records 
observers’ responses, and performs post-experiment data reduction. 

A chart of the sequence of computations performed during the presentation of a rotational trajectory 
stimulus in a typical psychophysical experiment is shown in figure 8. The prerun computations are 
shown in bold blocks and the relevant input parameters for each stimulus are shown in the circles. For 
each stimulus, the initial orientation bo with respect to the inertial system i, the direction and mag- 
nitude of the most direct path vector of rotation coo, and the total angle of rotation vo are specified. 

The orientation of bo in i is specified by the azimuth angle \|/, the elevation angle 0, and the roll 
angle <| ) (fig. 4). The direction of £)o is specified with respect to the bo-system by the azimuth angle y 
and the elevation angle 8 (fig. 5). The magnitude of coo is the specified rotational speed, coo- For each 
presentation the values of \|/,0,<J) and y,8 are picked at random and without replacement from a look-up 
table. 
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First, the quaternion Qi->b 0 is precomputed with vj/, 0, and <|>, using equation (17). The DCM D i->b 0 
is computed simply with equation (18). Next, following the geometry of figure 5, from y and 5, the 
components of fflo in the bo-system can be precomputed according to 


mo = ®o 


eye 5 
syc5 
-s5 


(42) 


Since fl>o = ffio coincides with the Euler axis for rotation from bo to bf and thus has the same direction, 
the fixed quaternion for rotation from bo to bf follows from equations (42) and (9): 


S.bo->bf 


coq/2 

cyc5soc(/2 

syc5sao/2 

-s8sao/2 


(43) 


and the corresponding quaternion matrix is computed with equation (14). The latter one, 

together with the quaternion of the initial orientation 9i->bo can use( ^ to P recom P ute the quaternion 
of the final orientation 


9i— >bf “ {MbQ-^bfJH.i^bfl 


(44) 


and the corresponding quaternion matrix Mj^ bf is computed with equation (14). 

Last, D b() _ > q 0 is precomputed with y and 5, using equation (29); Dq j _ xi2 with ((id, using equa- 
tion (37) where (J>d is picked at random from a table; and £^3 is defined with equation (41). The 
matrix Q s is computed using equation (21), with £0 S computed with equation (23), and coo is 
computed with equation (42). 

The on-line computations during stimulus presentation are aimed at computing the “disturbed” vec- 
tor of angular rotation ffi^(t) with equation (40), which constitutes the vector of rotation at which the 
object proceeds rotating at time t. This requires the computation of (Q b (t) = cooe(t) with equations (26) 
and (27), which is the undisturbed vector of rotation for the remaining most direct path between orienta- 
tion b and bf, and the computation of various DCMs for transfc rmations between the b-, i-, qo-, qr, 
q 2 -, and q 3 -systems. It also requires the computation of the “ar gular distance” from b to bf, a(t) with 
equation (11). With a(t), the angle of disturbance v(t) is commuted with equation (28). The disturbed 
vector of rotation is used to compute Q(t) with equation (21), v/hich in its turn, is used to compute the 
rotational path, described by the quaternion qj_> b (t) by solving the differential equation of equation 
( 20 ). 
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It can be shown that the averaged rotational deviation from the most direct path, (3, and the duration 
of the rotation, normalized with respect to the duration of the most direct path, (tf - to)coo/a<), depend on 
the parameters vo, e, and ao only, and are not affected by the initial orientation, the orientation of the 
Euler axis, or the angle of the plane of disturbance, (J>d. Hence these values can be precomputed and 
tabulated. For each set of vo, e, and ao, the rotational path is preexecuted twice: the first time to com- 
pute the time interval (tf - to) and the second time, using this interval, to compute the averaged devia- 
tion, p. 

At run time the relevant parameter for specifying the magnitude of distortion of the rotational path is 
P rather than_yo- Therefore, before starting the run, the value of vo for realizing the trajectory with 
the specified P is computed by linear interpolation from the tabulated values. In order to verify that the 
specified trajectory was executed, the actual rotational deviation P is computed on line, and at the end 
of the presentation of the rotational path, the averaged value is compared with the specified one. 


4. CONCLUSION 


This report has presented a systematic method for creating suboptimal rotational trajectories. Two 
parameters of the distortion metric, Vo and e, independently control the magnitude of the distortion and 
its “shape” (i.e., how the distortion is distributed over the time course of the trajectory). These sys- 
tematic disturbances can be introduced into a motion-control program to examine human observers’ sen- 
sitivity to rotational trajectory optimality. In addition, this algorithm can be used in other applications 
requiring systematic and reproducible disturbances of rotational trajectories. 
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Figure 2. The effect of the exponent e on the trajectory shape. 
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Figure 3. Method for computing the averaged deviation from the most direct path. 
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z q 0 

Figure 6. Definition of the qj-system. 


20 



Figure 7. Definition of the q2 and q 3 -systems. 
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Figure 8. Concluded. 
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